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Abstract 

55 . We study the role of interference in the process of quantum cloning. We show that in order to 

pg \ achieve better than classical cloning of a qubit no interference is needed. In particular, a large 

class of symmetric universal 1^ 2 qubit doners exists which achieve the optimal average fidelity 
5/6 for such machines, without using any interference. We also obtain optimal average fidelities for 
interference-free cloning in asymmetric situations, and discuss the relation of the quantum doners 

i-Q ' found to the Buzek-Hillery quantum doner. 
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I. INTRODUCTION 

The celebrated no-cloning theorem [l| states that unknown quantum states cannot be per- 
fectly copied, i.e. no quantum mechanical evolution exists which would transform a quantum 
state ip according to ip ^ 4' ^ 4' fo^' ^W unknown ip. This theorem is at the basis of the 
security of quantum key distribution systems, and has thus tremendous technological im- 
plications. From the theoretical point of view, it is a simple consequence of the linearity 
of quantum mechanical time evolution. A corresponding theorem for classical probability 
distributions can be proved just as easily: Propagating classical probability distributions 
linearly (e.g. phase space densities according to the Liouville equation) does not allow to 
transform any unknown probability distribution p according to p ^ p ^ p. So how could 
you download a perfect copy of this article? 

Classical and quantum information differ fundamentally by the fact that for classical 
systems we have (at least in principle) unrestricted access to the distributed quantities 
themselves (e.g. positions and momenta of particles in phase space), whereas for quantum 
systems we do not. In every-day language, we understand under classical copying measuring 
those quantities, and then preparing another system with the same values of those quantities. 
Indeed, a photocopy machine measures the position of bits of ink on the page, assumes that 
their momenta are zero (but could in principle also measure them), and then prepares an 
arbitrarily large number of other pages with the same positions of bits of inks with zero 
momentum. If we work with an ensemble of initial documents, differing slightly in the 
positions of bits of ink, such that it is worth talking about a probability distribution p, the 
latter will obviously be copied perfectly by the copy machine as well if it copies each member 
of the ensemble perfectly — but the transformation is manifestly non-linear. 

According to quantum mechanics, all information about a quantum system is encoded in 
the wavefunction ip, and ip never attributes a sharp value to both positions and momenta of 
the particles involved. More importantly, the apparent failure of "local hidden variable" de- 
scriptions of Quantum mechanical correlations confronted with experiments testing Bell- type 
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inequahties [2, y, m 15[ suggests that in a quantum mechanical system physical observables 
do not even have an objective value till the observable is measured. Quantum cloning aims 
therefore directly at reproducing the wavefunction rather than, say, the coordinates and 
momenta of individual particles. If we want to compare classical and quantum copying on 



equal footing, we should consider directly linear transformations (stochastic maps) of the 
initial probability distributions in the classical case, and this is what we are going to do 
below. 

Shortly after the publication of the no-cloning theorem, Buzek and Hillery (BH) showed 
that pretty good approximate cloning is possible [6]. They invented a 1 — > 2 quantum 
copying machine which takes a first qubit A in an unknown pure state \ip) = a\0) + (3\1) 
and a second qubit i? in a fixed known initial pure state |0) (an "empty page") as input and 
produces two mixed states p'^ = p'^ = ^\ip){ip\ + ^\4''^){4''^\, where \ip'^) = a*\l)—P*\0) is the 
state orthogonal to \ip). A single ancilla qubit is necessary to perform the transformation. 
The BH machine is manifestly symmetric, i.e. the fidelities Fj obtained as overlap of the 
final mixed states with the initial pure state of A, Fj = {iplp'Aip) , j = A,B, are the same. 
It is also universal, in the sense that it copies all initial states {ip) with the same fidelity 
Fj = 5/6. Later it was shown that 5/6 is indeed the optimal fidelity for symmetric universal 



1^2 doners of a qubit 







8j, and that this is the largest fidelity compatible with the 
non-signaling constraint imposed by special relativity [S]. 

For non-universal doners, the fidelities FjiO^cj)) = {ilj(9,(l))\p'j\ip{9,(f))) depend on the 
initial states 1^(0,0)), parametrized by polar and azimuthal angles 6 and (p on the Bloch 
sphere. It is then convenient to introduce average fidelities Fa = Fb, 

F=— del) deFj{e,(l)) sine. (1) 
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The best possible classical cloning of a qubit with symmetric average fidelities leads to 
Fa = Fb = 2/3 [lO(]. The cloning is classical in the sense introduced above, i.e. the doner 
acts with a classical stochastic map on the vector of probabilities corresponding to the initial 
state, while it has no access to the quantum coherences. We will provide a new and simple 
demonstration of this result below. The optimal classical fidelities are very easily matched by 
quantum cloning. A universal symmetric quantum doner with fidelities F^ = Fb = 2/3 can 
be constructed trivially by measuring the first qubit in an orthogonal basis chosen randomly 

le Bloch sphere, and preparing the second qubit in the state found 



and uniformly over t 

in the measurement Hl]. Another trivial symmetric universal cloning machine leaves the 
original qubit unperturbed, produces the new one in a random pure state, and swaps the 
two with probability 1/2. It has fidelities Fa = Fb = 3/4. 

As always, when a quantum protocol offers a higher performance than the best possi- 



ble classical scenario, it is natural to ask which quantum resource is at the origin of this 



difference. "Entanglement" and "Interference" are two obvious candidates [12]. In ll|] it 
was argued that quantum doners without coherent interaction between original qubit and 
copy should have fidelities limited by the one obtainable from the second "trivial" cloning 
strategy described above. Brufi and Macchiavello studied the entanglement at the output 
of a universal quantum doner (QC) and found that bipartite entanglement between two 
outputs vanishes for A^ — > M quantum doners (A^ identical input qubits, M output qubits) 
in the "classical" limit M ^ oo (A^ fixed) of infinitely many copies, and in fact as soon as 
M > N + 2. [13[. The fact that classical cloning discards the quantum coherences whereas 
quantum cloning need not, leads naturally to the hypothesis that interference should play 
an important role in quantum cloning. Put the other way round, it seems likely that a 
doner which uses no interference at all should not be able to outperform classical cloning. 
Surprisingly, we show in this article that the contrary is true: A large continuous class of 
quantum doners exists which use strictly zero interference, but which outperform classical 
cloning, and in fact give the maximal possible average fidelity 5/6 in the symmetric case. 

Evidently, in order to study this question, the concept of interference needs to be made 
precise. A quantitative measure of interference was introduced in [14] , and was used to study 
interference in quantum algorithms. This led to the hypothesis that an exponential amount 
of interference is needed in any unitary quantum algorithm which offers an exponential speed 
up over the corresponding classical algorithm. The hypothesis was supported by further 
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numerical evidence in a study of disturbed versions of Grover's and Shor's algorithms 

n 

In pj| we demonstrated that the interference in a randomly chosen quantum algorithm is 



with high probability close to the maximum possible value, whereas in [17|] we showed that 
interference plays at most a minor role in the transfer of quantum information through spin 
chains. 

Below we briefly review the interference measure. We will then formulate quantum doners 

n 

very generally in terms of a dynamical matrix p^ , and rewrite the interference measure using 
the dynamical matrix. The average fidelities are linear functions in the matrix elements, 
and it turns out that vanishing interference leads to linear constraints. Finding the optimal 
average fidelity becomes an instance of semi-definite programming, which we solve both for 
given Fa or Fb-, and for Fa = Fb- The latter case leads to one specific quantum doner, 
very similar to the BH quantum doner, and in particular with the same optimal average 



fidelities. We show, however, that starting from that machine a whole continuous class of 
other machines can be constructed with the same optimal average fidelities and vanishing 
interference. 

A final remark is in order concerning the applicability of our results to quantum broad- 
casting. Sometimes the term "quantum cloning" is restricted to machines which leave the 
outputs in a product state. "Quantum broadcasting" does not make any such restriction 
[l9|. In the following we will not make this distinction, i.e. we do not restrict ourselves to 
machines with factorizing outputs such that these machines might as well be called quantum 
broadcasters. 

II. INTERFERENCE MEASURE 

The main physical intuition behind the interference measure is that it should measure the 
coherence of a propagation, as well as equipartition: A classical stochastic map (i.e. a process 
which is not coherent at all), should not give rise to interference, whereas a unitary evolution 
might. However, a pure permutation of basis states would normally not be considered as 
creating interference, even if completely coherent. Clearly, initial states need to be "split" 
and superposed. The more states contribute with appreciable amplitude to a final state, 
the larger the interference. This is the "equipartition" property. In [14] we derived an 
interference measure by studying the change of final state probabilities as function of phases 
in the initial states. This led to an interference measure which, while probably not unique, 
does measure both coherence and equipartition. It maps any propagator P of a density 
matrix p, p' = Pp, to a real number between and A^ — 1, where A^ is the dimension of the 
Hilbert space on which p and p' act. Note that P is a super-operator, which, when specified 
in the computational basis {|i)}, propagates p with matrix elements pij = {i\p\j) according 
to 

Pij = 2_^ Pij.kiPki ■ (2) 

k^l j 1 

In terms of P, the interference X is given by [ij] 

^{P) = 2_^ \PiiM\ " /_^ \Pii,kk\ ■ (3) 

i,k,l i,k 

Clearly, if P only propagates probabilities (i.e. reduces to a classical stochastic map), Ph^m oc 
Skiik, I, where 6ki is the Kronecker-delta, X(P) vanishes. The squares in ([3]) assure the 



measurement of equipartition. This becomes most obvious for unitary propagation, p' = 
UpW , in which case I reduces to X = f A^ — '^- ^ \Uik\^ ) , where the double sum is nothing 
but the sum of inverse participation ratios of the column vectors of U. 

III. QUANTUM CLONING IN TERMS OF THE DYNAMICAL MATRIX 
A. Dynamical matrix 



ate 1— >■ 2 qubit quantum cloning very generally in terms of the dynamical 
120( 1 ■ D is related to P by a simple reshuffling of indices, D^j = Pikji, ^^^ 



We formu 
matrix D 

therefore contains all information about the propagation of the two-qubit system (original 
qubit and copy). In terms of D, eq.([2]) reads 



p'ij = ^D^^.Pki- (4) 



kJ 



The advantage of this is that if one considers ik and jl as single indices / and J (taking values 
0, . . . , 15 for ik = 00, 01, . . . 03, 10, 11, ... , 33, and similarly for jl), Dik = Djj, the matrix D 
acquires certain useful properties, inherited from the properties of p'. It can be shown that 
D is hermitian {Djj = Dj^) and block positive {{ip\D\ip) > Wl^p) = \u) ® \v), where \u) 
and |f ) are arbitrary pure single qubit states). Block-positivity implies immediately that all 
diagonal matrix elements of D must satisfy < Dij < 1. Furthermore, 

3 

Vd..=1Vj (5) 

i=0 



assures the correct normalization of p' [18|]. Given the linear nature of the propagation of 
the density matrix, eq.(jl]), the average over all initial states for the fidelities Fa and Fb can 
be performed once and for all, and leads to 



(6) 



Fa 


= -{2D 00 + Doo + D02 + 2Dm + Dio + D12 

6 V 00 22 02 10 32 12 




+ 1^20 + D22 + 2D22 + DSO + D32 + 2D32 

20 00 22 30 10 32/ 


Fb 


= -(2D00 + Doo + D02 + Dio + D12 +2D12 
5 V 00 12 02 10 00 12 



+2D20 +D20 +D22 +D30 +D32 +2D32) . (7) 

20 32 22 30 20 32 / 



Therefore, only 12 matrix elements of D determine the average fidelity of the A and B 
clones. Equations fl6|71) can be rewritten in the form 

Fa = A-D, Fb = B-D (8) 

with two hermitian matrices A and B easily extracted from eqs. (16|7|) . The multiplication 
in eq.Q is understood as the scalar product between two 16 x 16 matrices, i.e. 

15 

A-D= J2 -^iK-DiK. (9) 

I,K=Q 

As for the interference, eq.([3]) is rewritten as 

I{D) = J2\D^^\'-J2\^t\'- (1°) 

In order to have zero interference, it is then clear that we must have 

D^k =Oyk^l. (11) 

il 

In other words, when D is written as a matrix (with indices /, K introduced above), the 
off-diagonal matrix elements of the diagonal blocks of D must vanish. 

B. Classical propagation 

As discussed in the Introduction, we understand under classical cloning of a probability 
distribution the propagation of the probabilities with a stochastic classical map, 

p[, = J2D^^Pkk, (12) 

^ ^ ik 

k 

i.e. only the diagonal matrix elements Djj contribute. As we have furthermore eq.dS]), we 
obtain the average classical fidelities 

Fi = -(3 + Doo + Dio - D02 - D12) (13) 

g V 00 10 02 12 / 

fI, = - (3 + Doo + D20 - D02 - D22) . (14) 

g V 00 20 02 22 / 

Since < Doo + Dio < 1, and < D02 D12 < 1 (as well as < Doo + D20 < 1 and 

00 10 02 12 00 20 

< Do2,D22 < 1 concerning F^) we obtain immediate bounds for F^ and F^, 

l<FiX<l. (15) 



The upper bound reproduces the one found in [10| . Thus, unknown pure states of a single 
qubit cannot be cloned classically with average fidelity better than 2/3. Note that the 
upper bound is saturated by a classical identity map {Dik = 6ik), which copies perfectly the 

ik 

probabilities, but not the quantum coherences. 

IV. OPTIMIZATION 

The average fidelities Fa and Fb are linear functions and therefore convex in the matrix 
elements Dtk. There are 12 real independent matrix elements Dij on the diagonal, and six 
4x4 sub-blocks with non-vanishing complex matrix elements in the upper block triangle, 
such that D contains 204 real independent variables. If we are interested in a symmetric 
doner with F a = Fb-, or if we want to find the maximum or minimum value of Fa for 
given Fb (or vice versa), we add another linear constraint and effectively reduce the number 
of variables by 1. The optimization is over all dynamical matrices satisfying the above 
mentioned constraints (hermiticity, block positivity and partial traces equal one, as well as 
vanishing off-diagonal matrix elements in the diagonal blocks in order to have vanishing 
interference). The large number of independent variables in the optimization problem calls 
for a numerical solution. Fortunately, the problem falls into the class of convex optimization 
problems, as both the function to be optimized and the allowed domain for D is convex: 
it is very easy to show that the set of matrices which are block-positive is convex, and 
all other constraints are linear, i.e. do not change the convexity of the allowed domain. 
Furthermore, we note that if we impose the stronger constraint of positivity, the problem 
reduces to an instance of semi-definite programming (i.e. optimization of a convex function 
over a positive hermitian matrix, with eventual additional linear constraints). Routines for 
solving semi-definite programming problems are readily available. While one might worry 
that imposing positivity instead of block-positivity leads to a lack of generality, we show that 
for a symmetric doner with zero interference the maximum allowed average fidelity 5/6 is 
reached in the space of positive matrices, such that extending the space of matrices to those 
which are block-positive but not positive does not improve the result. In the following we 
will call T) the convex set of all positive dynamical matrices giving rise to zero interference. 

We used the "sedumi" routine of the YALMIP package under Matlab to solve the semi- 
definite programming problem, and found that for the symmetric case Fa = Fb the 
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dynamical matrix Z}°p* whose only non-vanishing matrix elements D'j^ are the diagonal 

Hiag- ('2iniililliiini^i'l the elements D°P* - D°P* - D°P* - D°P* - 
"-li^S 1,35 45 "J5 45 6' 4' 6' 4' 6' 4' 6' 4' ' 4' 3' 4/' eieuienis -L7o6 ~ -^^60 ~ -'^O 10 ~ -'^lOO ~ 

^4u = ^14 4 = ^su = ^14 8 = |, and Dlf = Dgf = D'^f^ = D°^l = i, leads to max- 
imum average fidelities Fa = Fu = 5/6. This is the maximum average fidelity possible 
for a symmetric universal doner 7|, |8|, |9(] . In particular, the maximum average fidelities are 
larger than the allowed classical value 2/3, and we have therefore demonstrated that better 
than classical quantum cloning is possible without using interference. Further below we will 
investigate this machine in more detail, and compare it in particular to the BH QC It is 
easy to check that the found QC is universal, i.e. all pure initial states are copied with the 
same fidelity Fa = Fb = 5/6. The minimal possible average fidelity of Fa for unrestricted 
Fb (and vice versa) is found to be 1/3. 

In the case of given Fa (or Fb), we show in FiglHthe maximum and minimum averaged 
fidelities Fb (or Fa) as a function of Fa (or Fb) for the full allowed range of the latter. We 
see that Fa can reach unity, if at the same time the other qubit is completely randomized 
{Fb = 1/2). The figure also shows that there is a whole continuous range of asymmetric QCs 
where both Fa and Fb are larger than the classical upper-bound2/3. The curves Fa {F'b) 
and Fb {Fa) enclose the domain af all possible QCs with positive dynamical matrix D. 
The domain is convex as Fa and Fb are linear functions of the Dik which form the convex 
domain V. 



V. MORE GENERAL INTERFERENCE-FREE SYMMETRIC QCS 

Semidefinite programming guarantees that a local optimum is also a global optimum. 
However, this does not exclude that other machines may exist with the same fidelity. Indeed, 
we shall show now that a high dimensional continuous class of symmetric QCs exists which 
all reach maximum average fidelity without making use of interference. 

First observe that the dynamical matrix Z}°p* has eigenvalues 1 (doubly degenerate), 
1/4 (8 times degenerate), and (6 times degenerate). Since Z}°p* is real and symmetric, 
its eigenvectors form a real orthogonal matrix O. The idea is then to add a hermitian 
perturbation W to Z}°p* in the subspace of eigenvectors of Z}°p* corresponding to the non- 
vanishing eigenvalues, and which leaves the partial traces unchanged. If the perturbation 
is small enough, the positivity of Z)°p* cannot be jeopardized immediately. Furthermore, a 
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FIG. 1: (Color online) The maximal and minimal average fidelity Fb of an interference free doner 
as function of F^ (black line) enclose the convex domain of all possible interference-free doners 
with positive dynamical matrix D. The maximum allowed classical fidelity for a symmetric doner 
is 2/3 (red lines). 

perturbation can be easily constructed such that it is orthogonal to both A and B in the 
sense of the scalar product IQ. Let wu be the matrix elements of VT (0 < /, J < 15) in 
the eigenbasis of D°^^. Transformed back to the computational basis, W = OWO^, we find 
A ■ W=B ■ W=5{woo + Wii)/12. Thus we have unmodified success probability Fa for 



Wu 



'Wqq. 



From the requirement of unchanged partial traces we get 

Wqo = Wii = 
W33 + W55+W77 + Wgg = 
^22+^44 + ^66 + ^^88 = 0. 



(16) 

(17) 
(18) 
(19) 



The only thing which remains to be imposed is that transformed back to the computa- 
tional basis the perturbation still has vanishing off-diagonal matrix elements on the di- 
agonal blocks, in order to keep the interference at zero. This requires wqi, wq2, W03, 
Wo4, Wo5, Woe, Wot, W14, W15, Wiq, wn, Wis, Wig, W23, W45, Wgj, Wsq to vanish. We 
are left with a large class of dynamical matrices Z}°p* which can be written as a 4 x 4 
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matrices of 4 x 4 complex sub-blocks D^'j with D^j 



Dfj and J, J G {0, ... , 3} 



where D^l^ = diag (|, i + ^99, 0, ^ + w^ss), ^ii* = diag (|, i + W77, |, i + w^ee), D 
diag (|, I + W55, 14 + ^44), D°f = diag(0,i + w^33, i'i + ^22) and 
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They depend on 64 real parameters, which, when taken small enough, leave iI)°P* positive. We 
have thus found a large continuous class of symmetric interference-free QCs with maximum 
average fidelity 5/6. It turns out that all these machines are also universal. This is a 
consequence of imposing zero interference, as is easily checked by relaxing this condition 
(see the discussion of vanishing matrix elements Wij after eq. 019p ). 

VI. COMPARISON WITH THE BUZEK HILLERY QC 

The comparison with the BH QC is not as straightforward as it may seem. The reason 
is that BH only specify what happens to a state \il)) ® |0), but the fate of \il)) ® |1) remains 
unspecified. In other words, only part of the dynamical matrix is defined, and in order to 
compare our interference-free QCs with the one by BH on the basis of the dynamical matrix, 
we would have to extend the definition of the BH QC to deal with the input state |1) for 
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the second qubit as well. If we choose a natural extension, 

\0)a\l)b\Q). -^ \0)a\l)b\Q2)x + {\0)a\0)b+\l)a\l)b)\Y2)cc 

\l)a\l)b\Q). ^ \l)a\0)b\Q3). + ma\0)b + \l)a\l)b) {Y^). , (23) 

where the subscript x denotes states of the doner, unitarity imposes the same constraints 
as in eq.(3.3) in [g], but extended to i = 0, . . . , 3. With an appropriate choice of overlaps 
between the states of the doner, one can show that the BH QC is a member of the continuous 
class of QCs with dynamical matrices Z}°p* found above. Alternatively, we can compare the 
QCs on the basis of the output reduced density matrices p'^ and p'^ obtained for mapping 
states lip) |0). Doing so, we find that all optimal symmetric interference-free QCs with 
dynamical matrices D°p* give the same single-qubit reduced matrix as the BH QC, 

^a p g 3 \a\ 



with p'^ = Tib [p'ab], p'b = Tr^ [p'ab] and pUb = I^°p*|^)(^| ® |0)(0|. 

VII. CONCLUSION 

We have found a large class of 1 ^ 2 qubit doners, which outperform classical symmetric 
cloning of qubits in terms of the average fidelities, but do not use any interference in the 
effective map of the two qubits. In particular we find QCs which reproduce the optimal 
average fidelities 5/6 of symmetric universal QCs without using interference. The well- 
known Buzek-Hillery QC turns out to be a member of this general class of interference- 
free QCs. How is it possible that interference-free quantum cloning outperforms classical 
quantum cloning? The answer is clear from the structure of the dynamical matrices D: 
Classical cloning corresponds to a diagonal D, quantum cloning to a potentially full matrix 
D. Interference-free quantum cloning is situated somewhere in between — only the off- 
diagonals of the diagonal blocks need to vanish. In other words, interference-free propagation 
of the density matrix does not map any coherences (i.e. off-diagonal elements of the density 
matrix) to probabilities, but may well map coherences to other coherences. The latter process 
influences the fidelities, and the corresponding matrix elements can be used to optimize the 
fidehties beyond classically possible values. It is intriguing to find that for optimal quantum 

12 



cloning (in the case of a 1 ^ 2 qubit doner) just this happens: only probabilities are mapped 
to probabilities, whereas coherences never modify the probabilities. It would be interesting 
to find out if this generalizes to higher dimensional or multiple-copy quantum cloning. It 
is also worthwhile noting that we have considered here the interference in the effective 
propagation of the two qubit system (original and copy). More general settings could be 
envisaged, for example including the copying machine as well. Additional information about 
the copying machine is needed then, however, and the results may depend on the dimension 
of the Hilbert space of the machine itself. 
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